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On dynamical gluon mass generation 
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Abstract. The effective gluon propagator constructed with the pinch technique is governed by a Schwinger- 
Dyson equation with special structure and gauge properties, that can be deduced from the correspondence 
with the background field method. Most importantly the non-perturbative gluon self-energy is transverse 
order-by-order in the dressed loop expansion, and separately for gluonic and ghost contributions, a property 
which allows for a meanigfull truncation. A linearized version of the truncated Schwinger-Dyson equation is 
derived, using a vertex that satisfies the required Ward identity and contains massless poles. The resulting 
integral equation, subject to a properly regularized constraint, is solved numerically, and the main features 
of the solutions are briefly discussed. 

PACS. 12.38.Lg Other nonperturbative calculations - 12. 38. Aw dynamics, confinement, etc 
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It is well-known that one of the main theoretical prob- 
lems when dealing with Schwinger-Dyson (SD) equations 
is that they are built out of unphysical off-shell Green's 
functions; thus, the extraction of reliable physical informa- 
tion depends crucially on delicate all-order cancellations, 
which may be inadvertently distorted in the process of 
the truncation. The truncation scheme based on the pinch 
technique (PT) [IJ[2] implements a drastic modification at 
the level of the building blocks of the SD series. The PT 
enables the construction of new, effective Green's func- 
tions endowed with very special properties; most impor- 
tantly, they are independent of the gauge-fixing parame- 
ter, and satisfy QED-like Ward identities (WI) instead of 
the usual Slavnov- Taylor identities. The upshot of this ap- 
proach would then be to trade the conventional SD series 
for another, written in terms of the new Green's functions, 
and then truncate this new series, by keeping only a few 
terms in a "dressed-loop" expansion, maintaining exact 
gauge-invariance. Of central importance in this context is 
the connection between the PT and the Background Field 
Method (BFM) , a special gauge-fixing procedure that pre- 
serves the symmetry of the action under ordinary gauge 
transformations with respect to the background (classi- 
cal) gauge field A a . As a result, the background n-point 
functions satisfy QED-like all-order Wis. The connection 
between PT and BFM, known to persist to all orders (last 
two articles in :2j), affirms that the (gauge- independent) 
PT effective n-point functions coincide with the (gauge- 
dependent) BFM n-point functions provided that the lat- 
ter are computed in the Feynman gauge. In this talk we 
report recent progress on the issue of gluon mass genera- 
tion in the PT-BFM scheme [3]. 



We first define some basic quantities. There are two 
gluon propagators appearing in this problem, A fll/ (q) and 
A^ v {q), denoting the background and quantum gluon prop- 
agator, respectively. Defining P fJ/V (q) = — , we 



have that A^ u (q), in the Feynman gauge is given by 



A^(q) 



P^(q)A(q 2 ) 



(i) 



The gluon self-energy, II ^(q), has the form il^ v {q) = 

P^{q)n(q 2 ), and A-\q 2 ) = q 2 + ifl(q 2 ). Exactly anal- 
ogous definitions relate A llv {q) with II^ v (q). 

As is widely known, in the conventional formalism the 
inclusion of ghosts is instrumental for the transversality 
of II^,(q), already at the level of the one- loop calculation. 
On the other hand, in the PT-BFM formalism, due to 
new Feynman rules for the vertices, the one-loop gluon 
and ghost contribution are individually transverse [5]. 

As has been shown in [3 , this crucial feature persists at 
the non-perturbative level, as a consequence of the simple 
Wis satisfied by the full vertices appearing in the diagrams 
of Fig.©, defining the BFM SD equation for A^(q) @]. 
Specifically, the gluonic and ghost sector are separately 
transverse, within each individual order in the dressed- 
loop expansion. 

Let us demonstrate this property for graphs (ai) and 
(a 2 ), given by 
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n ab (a)\ =- 



^AZ(kW b u' p X aAll,{k + q), 



[dk\r^: A a J{k) 



(2) 
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Fig. 1. The SD equation for the gluon propagator in the BFM. 
All external legs (ending with a vertical line) are background 
gluons, wavy lines with grey blobs denote full-quantum gluon 
propagators, dashed lines with grey blobs are full-ghost propa- 
gators, black dots are the BFM tree-level vertices, black blobs 
are the full conventional vertices, and white blobs denote full 
three or four-gluon vertices with one external background leg. 

where [dk] — d d k/(2n) d with d = A — e the dimension of 
space-time. By virtue of the BFM all-order WI 



(3) 



and using the tree-level r^ va p given in [5], we have 

i u n?ML = c A g 2 5 ab q , j[dk] A?(k) , 
9 v nfML = -CA9 2 S ab q, [{dk} A?(k) , (4) 



-n<$,{q)\ & ) = . The importance 



and thus, <f(fij*(q)\ 
of this transversality property in the context of SD equa- 
tion is that it allows for a meaningful first approximation: 
instead of the system of coupled equations involving gluon 
and ghost propagators, one may consider only the sub- 
set containing gluons, without compromising the crucial 
property of transversality. We will therefore study as the 
first non-trivial approximation for n^ v (q) the diagrams 
(ai) and (a2). Of course, we have no a-priori guarantee 
that this particular subset is numerically dominant. Ac- 
tually, as has been argued in a series of SD studies, in the 
context of the conventional Landau gauge it is the ghost 
sector that furnishes the leading contribution [5] Clearly, 
it is plausible that this characteristic feature may persist 
within the PT-BFM scheme as well, and we will explore 
this crucial issue in the near future. 

The equation given in ([2]) is not a genuine SD equation, 
in the sense that it does not involve the unknown quantity 
A on both sides. Substituting A — > A on the RHS of © 
(see discussion in [3]), we obtain 



[dk] rfA{k)w v 



i0 A(k + q) 



C A g 2 d 9liv [dk]A(k), 



(5) 



with = (2k + q)fig a p - 2q a g f _ L p + 2qpg m , and 



q v W va fi 



A- 1 {k + q)-A-\k) 



gap ■ 



(6) 



We can then linearize the resulting SD equation, by re- 
sorting to the Lehmann representation for the scalar part 
of the gluon propagator pQ 



A(q Z ) = d\ 



P(A 2 ) 



(7) 



(8) 



q 2 — A 2 + ie ' 

and setting on the first integral of the RHS of Eq. ([5]) 

where r^ a g must be such as to satisfy the tree-level WI 

Q V r^ a0 = [(k + qf - A 2 ] g aP - (k 2 - \ 2 )g a p . (9) 
We propose the following form for the vertex 



r, 



VOL ft 



r va p + Ci (2k + q) v + ^ [k 2 - (k + q) 2 



ga/3 



C3 



C2 

2<r 



[(k + q) 2 + k 2 ]y qp9 

vol Q_ctQvf3 )(io) 



which, due to the presence of the massless poles, allows 
the possibility of infrared finite solution. 

Due to the QED-like Wis satisfied by the PT Green's 
functions, A" 1 (q 2 ) absorbs all the RG-logs. Consequently, 
the product d(q 2 ) = g 2 A(q 2 ) forms a RG-invariant (/i- 
independent) quantity. Notice however that Eq.© does 
notencode the correct RG behavior: when written in terms 
of d(q 2 ) it is not manifestly g 2 -independent, as it should. 
In order to restore the correct RG behavior we use the 
simple prescription proposed in [1 , whereby we substitute 
every A(z) appearing on RHS of the SD by 



A(z) 



g 2 A(z) 
g 2 {z) 



= [l + bg 2 \n(z/n 2 )]A(z). 



(11) 



Then, setting b = ^Ca 
we finally obtain 



6 (ci+c 2 ) 

5 ' ' 



4+4ci+3c 2 
5 



d-\q 2 ) = q 2 ^K' + b 

- f1 2 / A 
+ ■yb dz z 
Jo 



9 2 /4 



dz \1 

Az 

1 K 



4z\ 1/2 d(z) 



T J g 2 (z) 

1/2 d(z) 



t{z) 



d-\o), 



K' 



and 



dz 1 + 7 ■ 



bo 



Az 
777 



(12) 

1/2 M 



d-i(0) = -^ d 4 k 



d(k 2 ) 
g 2 (k 2 ) 



(13) 
(14) 
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It is easy to see now that Eq. (fl2|) yields the correct UV 
behavior, i.e. d^ 1 (q 2 ) — b q 2 \n(q 2 / A 2 ). 

When solving (fT2]) we will be interested in solutions 
that are qualitatively of the general form p] 



d(q 2 ) = 



_2 



(<z 2 ) 



q 2 + m 2 {q 2 ) ' 



where 



—2 

9np 



(1 2 ) = 



b In 



q 2 + f(q 2 ,m 2 (q 2 )) 
A 2 



(15) 



(16) 



glsviq 2 ) represents a non-perturbative version of the RG- 
invariant effective charge of QCD: in the deep UV it goes 
over to g 2 (q 2 ), while in the deep IR it "freezes" [l][Z]j due 
to the presence of the function f(q 2 ,m 2 (q 2 )), whose form 
will be determined by fitting the numerical solution. The 
function m 2 (q 2 ) may be interpreted as a momentum de- 
pendent "mass" . In order to determine the asymptotic be- 

we 

1. 



havior that Eq. (fP2"]) predicts for m 2 (q 2 ) at large q 2 
replace Eq. (fT5"|) on both sides, set (1 — Az/q 2 ) 1 / 2 - 
obtaining self-consistency provided that 



m 2 (q 2 ) ~ ml ln" Q (q 2 / A 2 ) , with a = l+ 7 > . (17) 



The seagull-like contributions, defining d _1 (0) in 
are essential for obtaining IR finite solutions. However, the 
integral in (fT4|) should be properly regularized, in order to 
ensure the finiteness of such a mass term. Recalling that 
in dimensional regularization J[dk]/k 2 = 0, we rewrite the 
Eq.|l4|) (using (fTS])) as 



VlAk 2 ) 



[k 2 + m 2 (k 2 ))g 2 {k 2 ) k 2 
m 2 (k 2 ) 



k 2 [k 2 + m 2 (k 2 )} 

f{k 2 ,m 2 {k 2 )) 



-— \[dk]d{k z ) In ( 1 



The first integral converges provided that m 2 (k 2 ) falls 
asymptotically as ln~ a (fc 2 ), with a > 1, while the second 
requires that f(k 2 ,m 2 {k 2 )) should drop asymptotically at 
least as fast as ln _c (fc 2 ), with c > 0. Notice that perturba- 

tively e? -1 (0) vanishes, because m 2 (k 2 ) = to all orders, 
and, in that case, / = also. 

Solving numerically Eq. (fl~2|) . subject to the constraint 
of Eq.(fT4|). we obtain solutions shown in Fig.®; they can 
be fitted perfectly by means of a running coupling that 
freezes in the IR, shown in Fig. ([3]), and a running mass 
that vanishes in the UV [3]. er is treated as a free param- 
eter, whose values are fixed in such a way as to achieve 
compliance between Eas. p^l - dT^l) . 

Acknowledgments: This research was supported by Span- 
ish MEC under the grant FPA 2005-01678 and by Fundagao 
de Amparo a Pesquisa do Estado de Sao Paulo (FAPESP- 
Brazil) through the grant 05/04066-0. J. P. thanks the or- 
ganizers of QNP06 for their hospitality. 




Gluon Propagator 
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Fig. 2. Results for d(q 2 ), for different values for d _1 (0) (all 
in GeV 2 ), and the corresponding values for a. 




j-^g^j Fig. 3. The running charge, a(q 2 ), corresponding to the gluon 
propagator of Fig.©. 



References 



J. M. Cornwall, Phys. Rev. D 26, 1453 (1982). 
J. M. Cornwall and J. Papavassiliou, Phys. Rev. D 40, 3474 
(1989); J. Papavassiliou, Phys. Rev. D 41, 3179 (1990); 
D. Binosi and J. Papavassiliou, Phys. Rev. D 66, 111901 

(2002); J. Phys. G 30, 203 (2004).^ 

A. C. Aguilar and J. Papavassiliou, arXiv:hep-ph/0610040 
R. B. Sohn, Nucl. Phys. B 273, 468 (1986); A. Hadicke, 
JENA-N-88-19. 

L. F. Abbott, Nucl. Phys. B 185, 189 (1981). 
L. von Smekal, R. Alkofer and A. Hauck, Phys. Rev. Lett. 
79, 3591 (1997); R. Alkofer and L. von Smekal, Phys. Rept. 
353, 281 (2001); C. S. Fischer, J. Phys. G 32, R253 (2006); 
A. C. Aguilar, A. A. Natale and P. S. Rodrigues da Silva, 
Phys. Rev. Lett. 90, 152001 (2003); A. C. Aguilar, A. Mi- 
hara and A. A. Natale, Phys. Rev. D 65, 054011 (2002); 
Int. J. Mod. Phys. A 19 (2004) 249. 



